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INTRODUCTION 


Estimating the accuracy of stresses and deflections computed by the 
Finite Element Method (FEM) is a common problem that arises with each new 
application of the method. A computer code is now available for estimating 
the error in results from the STAGS finite element code (Ref. 1) for a shell 
unit consisting of a rectangular orthotropic plate. This memorandum 
describes the connection between the input data for the STAGS '’shell units” 
and the input data necessary to run the error analysis code (called EAC) . 
The STAGS code returns a set of nodal displacements and a discrete set of 
stress resultants. The continuous solution is defined by a set of 
generalized coordinates computed in EAC. The theory and the assumptions that 
determine the continuous solution are also outlined in this memorandum. 

A structure, such as a stiffened panel, is modeled in STAGS as a number 
of shell units and/or finite elements. The current version of EAC performs 
error analysis on results for a rectangular plate unit from STAGS. The 
theory in the error analysis requires input data for the plate geometry, the 
material properties, and a set of boundary conditions. These input data 
come from the STAGS code. The specific link between the input and output 
data from STAGS and the input data for EAC is POSTP, the postprocessor 
program in the STAGS processors . A modified version of POSTP is described in 
Section 2 of this memorandum, following a brief review of the general 
procedure for executing EAC as part of the STAGS processors in Section 1. 

Section 3 presents the theory for an orthotropic, rectangular plate with 
initial geometric imperfections. The analysis for the continuous solution is 
outlined, and the notation for the output data from EAC is defined. The 
error analysis is for a postbuckled plate. The postbuckled finite element 
solution in STAGS can be computed using the new TP option in STAGS (Ref. 2). 
Input for running the TP option in STAGS is summarized in Section 4. 
Finally, input and output for a sample problem are contained in Appendix B. 
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SECTION 1: STAGS DATA 


Most of the input data for EAC is STAGS input or output data. The STAGS 
processors consist of STAGS1, STAGS2 , POSTP, and STAPL. The input data for 
the plate geometry are part of the STAGS1 input. The material properties and 
lay-up for an orthotropic plate are also part of the usual STAGS1 input data. 
The error analysis in EAC uses STAGS2 output for edge displacements and 
rotations as boundary conditions. A complete description of STAGS input is 
found in the STAGS Users Manual, Ref. (1). Only input directly connected to 
the error analysis is reviewed here, the input data for the finite element 
results in STAGS are unchanged from the current version of STAGS . 

1.1 Plate Geometry 

In STAGS a structure, such as a stiffened panel, is modeled as a set of 
"shell units" and/or a set of finite elements. An execution of EAC is 
concerned with a single shell unit, a rectangular plate. The geometry of the 
plate is defined as STAGS1 input. The input follows the format starting with 
ISHELL=2 on the M-l record. The length and width of the plate are deduced 
from data on the M-2A record. The logic in EAC assumes that a standard 

rectangular grid of nodal points is defined for the plate in STAGS . The 
number of rows in the rectangular grid for shell unit I is NRWS(I) on record 
F-l, and the corresponding number of columns is NCLS(I). There are several 
options in STAGS for the definition of initial imperfections over a flat 
plate. Any of these options can be used in the FEM analysis, however, the 

logic in EAC assumes a continuous function for the imperfections. Therefore, 
input for the definition of imperfections is required both in the STAGS1 
input data and for EAC (see instructions for the EAC-6 record in Section 2), 

1.2 Stiffness Properties 

The theory in Section 3 of this memorandum defines stiffness 

coefficients A^., A 12 , A^, k ^ , D n> D^, D^, and D &6 for an orthotropic 

plate. These stiffnesses are computed by STAGS1 from the input in the K 
records and IWALL in the M-5 record. 
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1.3 Boundary Conditions 


The boundary conditions for the error analysis in EAC for a plate shell 
unit are displacements and rotations on the edges of the plate computed in 
STAGS2. The modification of the postprocessor POST? to manage these data in 
EAC is discussed in Section 2. The STAGS1 input data for boundary conditions 
and loads for the structure are unchanged; the specific plate unit considered 
in EAC has the usual input for boundary conditions, loads, and connection of 
boundary lines. 

After completing the input data for STAGS1 for the shell structure 
containing the plate unit, STAGS1 and STAGS2 are executed to compute a finite 
element solution for the shell structure. The error analysis in EAC is then 
executed as a postprocessor of STAGS2 data. The data that must be extracted 
for the plate unit in EAC are the geometry, the stiffnesses, and a set of 
boundary conditions. The details of the postprocessing to obtain input for 
EAC are contained in the next section. 
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SECTION 2: POSTPROCESSING STAGS DATA 


The error analysis code EAC requires a few new and some modified input 
data files and procedures. Detailed instructions are provided in this 
section to help the user prepare the various input files. Example procedure 
files are also provided for the execution of the STAGS postprocessors on the 
Cyber and MicroVAX machines at NASA Langley Research Center [LARC] . 

As mentioned earlier in Section 1, EAC uses as input, selected data from 
the input to and output from a STAGSC-1 analysis. This input file, called 
TAPE43, is prepared by using NPOSTPX, a modified version of the STAGS 
postprocessor POSTP. POSTP is used following a previous STAGSC-1 execution 
to enable the user to print quantities not selected for output in the 
original execution. NPOSTPX is POSTP modified to output selected data to a 
separate file, TAPE43. NPOSTPX requires a separate input file to specify the 
particular plate unit in the structure being analyzed, and the load step at 
which the data are to be written to TAPE43. 

In its present version, EAC can handle only one plate unit, at one load 
step of the STAGS execution, using the non-corotation option. The FEM 
solution from STAGS is assumed to be a valid zeroth approximation for 
starting the iterative solution of the shell equations. This solution forms 
the first step of execution of EAC, during the first phase of which the 
resulting general solution in the form of Fourier coefficients and constants 
of integration (generalized coordinates) are written to a local file called 
TAPE45. During the second phase, these data are read from TAPE45 and are 
used to generate results at specified points on the plate, or over a user 
defined grid. There exists an option in EAC which permits the user to reuse 
data on TAPE45 saved from a previous EAC execution to generate output data 
over a different output grid or at selected points on the plate unit. 
Section 3 contains a summary of the numerical analysis in EAC that results in 
the generalized coordinates for the continuous solution. Results are written 
to two local files - TAPE6 and TAPE46 . 

TAPE46 is the input data file for EACPL , a contour plotting post- 
processor for EAC. TAPE46 contains all the functional data for plotting two 



dimensional contour plots. Currently EACPL uses the USCONTOUR plotting 
package on the Cyber and DI3000 contour package on the MicroVAX. EACPL 
requires an additional input file called TAPE47 to provide plotting details 
such as the contour limits, intervals, etc. Examples for the MicroVAX are 
also provided in subsections titled WORKING ON THE MicroVAX. 

2.1 INPUT 

2.1.1 FEM INPUT 

EAC requires as input, a local file called TAPE43 . This is a file of 

selected discrete results from STAGS to be analyzed. TAPE43 is obtained by 

executing a modified version of POSTP after a previous STAGSC-1 execution. 
The input required for this new version of POSTP (called NPOSTPX) is 
described in the next sub-section. A sample procedure file for the 

generation of a TAPE43 is included in Section 2. STAGS results selected for 
TAPE43 are displacements, rotations, stress and moment resultants, strains, 
changes in curvatures. As with POSTP, the NPOSTPX execution requires a 
Solution Data file (SOD file) generated by a previous STAGS execution. On 
the Cyber, this file is called TAPE22 and must be available as a local file 
named TAPE20. A review of the POSTP input data cards as applicable to 

NPOSTPX is given next. 

Note: The equivalent of NPOSTPX on the VAX is POSTP.EXE 

2.1.2 INPUT TO NPOSTPX 

TAPE43 , the input data file for EAC, is created by executing NPOSTPX. 
NPOSTPX is a modified version of POSTP and was designed to create a TAPE43 as 
well as to select data created and saved by a previous STAGS execution for 
printing or plotting. It uses as input, a subset of the POSTP input. Just 
as with POSTP, NPOSTPX can be repeatedly executed with the same STAGS SOD 
data, each time with a different selection of output parameters. This would 
be useful in the preparation of TAPE43s for a series of EAC executions for 
different plate units in the same structure. POSTP users are directed to the 
STAGS Users Manual, Ref. 1, for input instructions. All input records to 
NPOSTPX are in STAGS free field format and reviewed below. 
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Note 1: EAC can only be used for one plate unit and one load step at a 
time. 

Note 2: The co-rotation option in STAGS is incompatible with EAC, because 
edge rotations from the co-rotational option are not defined with 
respect to the local fixed coordinate system. 

Note 3: When used solely for the creation of a TAPE43, the user needs to 
specify only the first four records PP-1 to PP-4, with the first 
four variables in the PP-4 record selected as unity for the 
plate unit under consideration, keeping all others zero. 

PP-1 Case Title Record 

Variable Description 

COMENT Arbitrary 80 character alphanumeric field 

to be printed as case title. 

PP-2 Number of Load Steps 

Variable Description 

NSTEPS Number of load steps for which 

postrocessing is desired. Currently, 

NSTEPS*! . 


PP-3 Load Step Sequence 

Variable Description 

ISTEP Integer load step in STAGS to be used for 

error analysis. 

The next three records must be repeated for each shell unit in the 
structure. However, data to be printed must be selected only for the shell 
unit for which the TAPE43 is required. 


PP-4 Output Control Parameters 


Variable 

PRD 

I PR R 


Description 

0 - Do not print displacements 

1 - Print displacements for the load step. 

0 - Do not print stress resultants 

1 - Print stress resultants 
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IPRE 

IPRS 

IPRP 

IPRF 

NS ELD 

NSELS 

IPRSDP 

IPRSTR 

ISL 


ISS 

ISD 


0 - Do not print strains 

1 - Print strains and changes of curvature 

0 - Do not print stresses 

1 - Print stresses 
0 - for NPOSTPX. 

0 - Do not print point or nodal forces 

1 - Print nodal forces 

The number of records defining selective displacements 
The number of records defining selected stresses. 

0 - Print selected displacements at the load step 
0 - Print selected stresses at the load step 

0 - Stresses, strains and resultants are 

computed at element centroids. 

1 - Stresses, strains and resultants are 

computed at element integration points. 

0 - This entry is zero always . 

0 - Print stresses or strain components in the 

directions of the shell. 

1 - Print stress components in both the shell 

and the material coordinate directions. 


NSELD >0, go to PP-5 

NSELD - 0, NSELS >0, go to PP-6. 


PP-5 Output Control Record 2 

This record defines a number of nodes at which displacements are to be 
printed at each load or time step. If IROWD is set equal to zero the 
displacements are printed for each row, that is, the entire column. If ICOLD 
is set equal to zero the displacements are printed for each column, that is, 
the entire row. The record is repeated NSELD times (PP-4) . 

Variable Description 

IROWD Row number at which displacements are to 

be printed. 

ICOLD Column number at which displacements are to 

be printed. 

NSELS >0, go to PP-6. 

NSELS « 0, follow instructions at end of PP-6. 
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PP-6 Output Control Record 3 

This record defines a number of locations at which stresses will be 
printed at each load or time step. Stress output will be consistent with the 
description for ISO on the PP-4 record or in a user-written subroutine. The 
record is repeated NSELS times (R-l record of STAGS input) . 

Variable Description 

IROWS Row number identifying the element for 

stress output. 

ICOLS Column number identifying the element 

for stress output. 

2.1.3 ADDITIONAL INPUT FILE 

In addition to TAPE43 , a separate user prepared input file called TAPE44 
is required to control the execution of EAC and its output. The input data 
consist of a sequence of free field format input records, one record 
corresponding to each FORTRAN input list. A free field format input record 
contains a number of numerical data fields, each separated by one or more 
blanks . 

EAC-1 Case Title 

The Case Title is read on one line that may contain any character data, 
not exceeding 80 characters in length. 

EAC -2 Execution Type Selection Record 

Quite often, users may wish to rerun EAC for the same plate unit but 
with output generated for a different set of points or for distinct points on 
the example shell unit. There is an option which permits repeated execution 
of EAC for the same plate using a TAPE45 file generated by a previous 
execution . 

Variable Description 

IRERUN 0 - Begin new case. Read EAC-2A,B 

1 - Repeat with existing TAPE45 . Go to EAC-3 
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EAC - 2 A Orthotropic stiffness coefficient record 1 

The user must input the stiffness coefficients for the corresponding 
plate unit. These coefficients are available in the STAGS1 output for each 
component unit of the shell structure, and the user should select the 
appropriate? data. The required membrane stiffness coefficients are 
Variables: All, A12, A22, A66 

EAC-2B Orthotropic stiffness coefficient record 2 

The required bending stiffness coefficients are 
Variables: Dll, D12, D22, D66 


EAC -3 Solution Summary Record 

The discrete data input to EAC via TAPE43 are stored in arrays whose 
dimensions are based on the STAGS grid. The displacement data arrays 
correspond to the nodes of this grid, while the stress and moment resultant 
data arrays correspond to the mid-points of the rectangular elements. The 
dimensions of these data arrays are limited by the core memory limitations on 
the host computer. These limits have not been established as yet. At 
present, the dimensions of the input arrays are limited to a maximum of 50 
rows and 50 columns. 

EAC can generate results either at a point (IPT > 0) or over a grid (IPT 
* 0) . In the latter case, the dimensions of the output: arrays, consisting of 
displacements, normal stress and moment resultants, transverse shears and an 
error estimate, correspond to the number of nodes in the grid. Here again, 
the dimensions are limited to a maximum of 50 rows and 50 columns. If 
necessary, the user can access a copy of the EAC source code and increase 
these limits. Instructions on how to do this are provided in Section 2. 


Variable 

IPT 

MT 

NT 

MP 


NQ 


Description 

0 - generate output results over user grid 
>0 - generate results at IPT points 
Number of terms in X for product solution 
of in-plane homogeneous equations. (See 3.6) 
Number of terms in Y for product solution 
of in-plane homogeneous equations. (See 3.6) 
Number of terms in X for product solution 
of transverse homogeneous equation. (See 3.5) 
Number of terms in Y for product solution 
of transverse homogeneous equation. (See 3.5) 
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If IPT«0 then go to EAC-4, 

If IPT>0 then read IPT EAC-3A cards, go to EAC-5. 

Note: See section 3.5 for limits on MT , NT, MP , NQ . 

EAC-3A Point Coordinate Data Record 

If IPT > 0 then provide IPT records with the X and Y coordinates of points 
for output data. There is a limit of 10 records, i.e., IPT < 10. 


Variable 

X 

Y 


Description 
X coordinate of point 
Y coordinate of point 


EAC - 4 Output Grid Definition Record 
Variable Description 

NXS Number of grid points in X-direction 

NY'S Number of grid points in Y-direction 

EAC-5 Output Summary Record 
Rudimentary output control is provided. 

Variable Description 

IPRT 0 - No output printed. 

1 - Output printed on TAPE6 

IGENC 0 - Generalized coordinates not printed 

1 - Print generalized coordinates 


EAC -6 Initial Imperfection Record 

The presence of initial imperfections are accounted for in EAC by an 
imperfection subroutine WIMPD . The imperfection displacements are of the 
form 

W = | =1 IMN(J) SIN M(J)ttX/A SIN N(J)ttY/B 

where IMN, M and N are limited to 24 each. 

Variable Description 

IWIMP 0 - No imperfections 

1 - Imperfection function parameters 

specified in records EAC-6A to EAC-6D 
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EAC-6A Imperfection Number Record 


Variable Description 

K Number of imperfections 

EAC-6B Imperfection Coefficient Record- 

Variable Description 

IMN K imperfection coefficients 

EAC-6C M(J) Data Record 

Variable Description 

M(J) K M(J) data entries 

EAC-6D U(J) Data Record 

Variable Description 

N(J) K N(J) data entries 

If IWIMP=0, the input data file is complete. 

If IWIMP-1, then read records EAC-6A to EAC-6D. 

See Appendix A for nomenclature and Section 3 for theory. 

2.1.4 INPUT TO PLOTTING POSTPROCESSOR EACPL 

EACPL is the contour plotting postprocessor for EAC. The contour 
information is generated by EAC and stored on TAPE46 . EACPL uses these data 
to plot two dimensional contour plots of the input and output of EAC. For 
general runs, the plots generated are of 

1. STAGS OUTPUT - Displacements, stress and moment resultants 

2. EAC OUTPUT - Displacements, stress and moment resultants, 
transverse shear, and error in the transverse equilibrium 
equation . 

For point analysis runs (IPT>0 on EAC-3 card above), no EAC output plot data 
are generated, and consequently these plots will not be generated. EACPL 
needs an additional input data file called TAPE47 . 
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EACPL-1 Grid Size Data Record 


Variable 

NRWS 


NCLS 


NXS 


nys 


PSCALE 


ICON 


Description 

Number of nodal points in the X-direction 
for the STAGS grid. 

(repeated from F-l card in STAGS1 data) 

Number of nodal points in the Y-direction 
for the STAGS grid. 

(repeated from F-l card in STAGS1 data) 

Number of nodal points in the X-direction 
for the user (output) grid. 

(repeated from EAC-4 card above) 

Number of nodal points in the Y-direction 
for the user (output) grid. 

(repeated from EAC-4 card above) 

Scale factor (based on plate dimensions) . 

A value greater than unity enlarges the 
plot and a value smaller than unity will 
scale down the plot. 

0 - system will select the contour details 

based on plot data. 

1 - contour details given by user on 12 EACPL-2 cards. 


EACPL-2 Contour Limit Detail Record 

The contour limit details are given for the 12 parameters plotted. These are 
in order: U, V, W, NX, NY, NXY, MX, MY, MXY, QX, QY, ERR. 


Variable 

SLO 

SHI 

SING 


Description 

lover limit on the function to be plotted, 
upper limit on the function to be plotted, 
step size, such that No. of Contours < 40. 


The plotting program will use the same limits for corresponding STAGS and EAC 
plots . 
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2.2 PROCEDURES 


The STAGSC-1 code is currently available on the Cyber machines for users 
at NASA Langley Research Center. A new experimental version of the STAGS 
code is currently installed on a MicroVAX in Building 1148. Given below are 
sample procedure files and instructions for the execution of STAGS and EAC on 
both these machines. Users are cautioned that the procedure files given 
below are not complete and must be treated as examples only. Users on other 
systems must consult their systems manager to interpret these instructions. 
The RESOURCE identification of the MicroVAX system on MICOM is 85. 

2.2.1 WORKING ON THE CYBER 

A complete execution of the STAGSC-1 system would typically involve 
executing three separate but independent modules; namely STAGS1, STAGS2 and 
STAPL. STAGS1 reads all input data describing the problem, preprocesses this 
input, writes local files containing data to be used by STAGS2, the main 
computational module, and STAPL, the plotting postprocessor . STAGS2 does not 
read any user -prepared input data, so a separate record in the input data 
deck is not included for STAGS2. STAPL, however, does require user-prepared 
input data describing plots to be generated and, if executed during a run, 
requires a record in the input data deck. The user may opt to bypass STAPL 
and obtain the plots using EACPL, the plotting postprocessor for EAC. EACPL 
generates plots of the displacements, stress and moment resultants for both 
the STAGS discrete output and the EAC continuous approximations. The 
advantage of using EACPL rather than STAPL is that the plots for the discrete 
variables and the continuous variables will both agree in scale and contour 
intervals for ease of direct comparison. In addition, transverse shears and 
the error in the transverse equilibrium equation computed by EAC are plotted. 

The fourth module, which is either POSTP or NPOSTPX, can be used 
following a previous STAGSC-1 execution. POSTP enables the user to print 
quantities not selected for output in the original execution, while NPOSTPX 
is required to be executed to create a TAPE43 file as input to EAC. In 
addition to the creation of TAPE43 , NPOSTPX enables the user to print output 
data for the plate unit and the load step under consideration. A new set of 
procedure files are provided to access copies of the absolute binary files 
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of STAGS 1 , STAGS 2 , STAPL, and POSTP/NPOSTPX or to create new absolute binary 
files incorporating any user-written subroutines that may be required. The 
I master procedure file, STAGS, is accessed from the user's control card record 

as follows: 

GET , STAGS/UN-STAGS , ST-CPF . 

This file contains seven procedures; four associated with standard 
STAGSC-1 executions and three associated with STAGSC-1 executions using user 
written subroutines. The four procedures associated with standard execution 
are LITTLE, BIG, POST, and TP. Both LITTLE and BIG set up local copies of 
the absolute binary files for STAGS1, STAGS2 and STAPL. LITTLE should be 
used for models having less than 3,000 degrees of freedom. To access one of 
these procedures requires either of the following control cards: 

BEGIN(LITTLE, STAGS) 
or 

BEGIN (BIG, STAGS) 

The procedure POST sets up local copies of the absolute binary files for 
STAGS 1 , POSTP. and STAPL. POSTP processing is performed using a local file 
named TAPE20 (TAPE22 previously) and creates a new local file TAPE22. STAPL 
postprocessing is performed using the local file TAPE22. POST is accessed as 
follows : 

BEGIN(POST, STAGS) 

If user-written subroutines are included, it is necessary to reload the 
i STAGSC-1 processors using one of the following procedures: RLBIG or RLPOST. 

To access these procedures requires the following control cards: 

BEG IN (RLBIG , STAGS , UWS-lfn) 
or 

BEGIN(RLPOST , STAGS , UWS-lfn) 

where the selection of which procedure to use is determined by the user's 
application and lfn is the name of a local file containing the FORTRAN-77 
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source statements of the user-written subroutines. If these reside on a 
permanent file, the user must GET or ATTACH this file and specify the local 
name lfn. The RL. . . procedures will compile Ifn and generate new local 
absolute binary files for STAGS1, STAGS2 and POSTP. 

In order to use POSTP/NPOSTPX, TAPE22 from a previous STAGS execution 
must be available as a local file named TAPE20, In addition, the user must 
have prepared and executed STAGS1 with the same model that produced the 
original TAPE22. Model data must be unaltered, but the print control flags 
and IPOSTO, IP0ST1, and IP0ST2 may be changed as explained in the user's 
manual. It is not necessary to execute all modules in a single run. 

In order to use EAC , TAPE43 generated by postprocessing with NPOSTPX 
must be available as a permanent file. Other data files needed are: 


TAPE44 

user input to control execution of EAC 
described above under the heading 
ADDITIONAL INPUT FILE. 

TAPE45 

file containing the solution in generalized 
coordinates, created by a prior EAC run. 

TAPE46 

input file for EACPL, generated by EAC 

if I RERUN-1 on EAC-2 card. (See ADDITIONAL INPUT) 

TAPE47 

input file to control execution of EACPL. 
(described as INPUT TO EACPL) 
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2. 2. 1.1 STANDARD BIG STAGSC-1 ANALYSIS 

JOB. . . 

USER. . . 

CHARGE. . . 

DELIVER. . . 

FILESET , IA-CPF , DA-LPF . 

GET , STAGS/UN-STAGS , ST-CPF . 

BEGIN(BIG, STAGS) 

STAGS 1. 

STAGS 2. 

REPLACE (TAPE22-=pfn) (* pfn is the permanent file name of TAPE22 *) 
STAPL. 

PLOT . VARIAN , FF 
E-O-R 

STAGS1 analysis data 
E-O-R 

STAPL plotting data 
E-O-F 
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2. 2. 1.2 STANDARD POST- PROCESSING 


This procedure file is given for reference only. The procedure file to 
execute NPOSTPX to generate TAPE43 for EAC is given next. 


JOB. . . 

USER. . . 

CHARGE . . . 

DELIVER. . . 

FILESET . IA-CPF , DA-LPF . 

GET , STAGS/UN-STAGS , ST-CPF . 
GET , TAPE20~pfn . 

BEGIN(POST, STAGS) 

STAGS 1. 

POSTP. 

STAPL. 

PLOT . VARIAN , FF 
E-O-R 

STAGS1 data 
E-O-R 

POSTP data 
E-O-R 

STAPL data 
E-O-F 
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2. 2. 1.3 POST-PROCESSING FOR GENERATION OF TAPE43 


This is a skeleton procedure to execute NPOSTPX on the Cyber. It 
requires a TAPE20 file which was previously saved as TAPE22 from a STAGS-C1 
execution. TAPE43 is required to analyse the STAGS results. 

JOB. . . 

USER. . . 

CHARGE. . . 

DELIVER. . . 

FI LESET, IA-CPF , DA-LPF . 

GET , STAGS /UN- STAGS , ST-CPF . 

GET , NPOSTPX/UN— STAGS , ST-CPF . 

GET ,TAPE20-pfnl . (* Where pfnl is a previously saved TAPE22 *) 

BEGIN( POST, STAGS) 

STAGS1 . 

NPOSTPX. 

REPLACE , TAPE43-pfn2 . 

E-O-R 

STAGS1 analysis data 
E-O-R 

NPOSTPX data 
E-O-F 
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2. 2. 1.4 STANDARD EAC EXECUTION 

JOB. . . 

USER. . . 

CHARGE. . . 

DELIVER. . . 

FILESET , IA-CPF , DA-LPF . 

GET , EAC , EACPL/UN-STAGS , ST-CPF . 

GET , TAPE43~pfnl . 

GET , TAPE44-pfn2 . 

GET,TAPE45=pfn3. (* ONLY FOR A RERUN *) 
GET , TAPE47-pfn4 . 

EAC , TAPE43 , TAPES . 

REWIND, TAPE6 , TAPE46 . 

REPLACE , TAPE6-= 

REPLACE, TAPE46- 

COPYSBF, TAPE43 . (* PRINT DATA FILE *) 

COPYSBF , TAPE6 . (* PRINT OUTPUT FILE *) 

EACPL. 

PLOT . VARIAN , FF 
E-O-F 
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2.2. 1.5 STANDARD BIG STAGSC-1 ANALYSIS 


If user-written subroutines are included then the following control 
cards will be required. 

JOB. . . 

USER. . . 

CHARGE . . . 

DELIVER. . . 

GET , STAGS/UN-STAGS , ST-CPF . 

BEGIN(RLBIG , STAGS , UWS-INPUT) 

STAGS 1. 

STAGS 2. 

STAPL. 

E-O-R 

User-written subroutines 
E-O-R 

STAGS 1 data 
E-O-R 

STAPL data 
E-O-F 


The user-written subroutines could be in the file 'pfn' , in which case 

use : 

BEGIN (RLBIG , STAGS , UWS-pfn) 
and delete the statements: 

E-O-R 

user-written subroutines 


I 
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2.2. 1.6 MODIFYING EAC FOR LARGE GRIDS 


At the present time, the input and output array dimensions are limited 

to a maximum of 50 rows and 50 columns each. To handle larger arrays, the 

user can access and modify a copy of the EAC and EACPL source codes. These 
are accessed by using the control cards: 

GET , EACNOW/UN-STAGS , ST-CPF . 

GET , EAC PLOT/UN-STAGS , ST-CPF . 

To modify the limits on the INPUT (STAGS) arrays (determined by the F-2 
card input in STAGS1) , edit the source file and modify the values of SNX and 
SNY in the PARAMETER statements. SNX and SNY correspond to NRWS and NCLS 
respectively. 

To modify the limits on the OUTPUT (USER) arrays, modify the values of 

NXU and NYU in the PARAMETER statements. For either the STAGS output arrays 

or the EAC output arrays, edit and modify the limits in EACPLOT. To compile 
and load, use the following procedure file: 

JOB. . . 

USER. . . 

CHARGE. . . 

DELIVER. . . 

FILES ET , IA-CPF , DA-LPF . 

GET, EACNOW, EACPLOT. 

GET , PL0T3D , GRAFL IB/UN- STAGS , ST-CPF. 

FTN5 , I-EACNOW , L-MINE , LO-S , DB-PMD . 

ATTACH , FTN5MLB /UN-LIBRARY . 

LDSET , LIB-FTN5MLB , MAP-S BEX/MINE . 

LOAD , LGO 
NOGO , EAC 
RETURN , LGO . 

FTN5 , I-EACPLOT , L-MINE , LO-S , DB-PMD . 

ATTACH, FTN5MLB , LRCGOSF/UN-LIBRARY . 

LDSET , LIB-FTN5MLB/PLOT3D/GRAFLIB/LRCGOSF , MAP-SBEX/MINE . 

LOAD , LGO . 
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NOGO , EACPL . 

REPLACE, EAC, EACPL. 

ENQUIRE. 

EAC and EACPL are the new executables. 

NOTE : Cyber has limited resident core memory, which restricts the field 

length available to the 'user. What this means is that there is a finite 
upper limit to which the grid size can be increased. This limit has not been 
determined. 
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2.2.2 WORKING ON THE MicroVAX 


An experimental version of the STAGSC-1 computer code was installed on a 
MicroVAX in Building 1148 during May 1987. Currently, this version of STAGS 
does not have the STAPL postprocessor installed. POSTP is available and the 
VAX version will write a TAPE43 file whether the user wants to or not. 

A copy of the error analysis code (EAC) has also been installed on the 
MicroVAX. System symbols have been defined to run STAGS and EAC on the 
MicroVAX and detailed instructions are given in the following pages. 

2.2.3 EXECUTING STAGS ON THE VAX 

The user is reminded that this is an experimental version of the 
STAGSC-1 code and its capabilities and reliability are yet to be documented. 
In this version, an additional input file is required for the execution of 
STAGS2 . To run STAGS, the user enters the command : 

RUNSTAGS <CR> 

Entering this symbol initiates a command file that will submit a STAGSC-1 
analysis as a batch job. The user must use the following naming convention 
for the input files and resulting output and restart files. 


data . inp 

defines the STAGS1 input data 

data.bin 

defines the STAGS2 input data 

data . inpp 

defines the POSTP input data 
(NOTE: different from Cyber!) 

data. sod 

defines the solution data file 

data . rs t 

defines the previous restart file 

data. outl 

defines the STAGS1 run output 

data . out2 

defines the STAGS2 output 

data . out3 

defines the RESTART output 
(different versions for restarts) 

data . outpp 

defines the POSTP run output 

data. t43 

defines the TAPE43 created by POSTP 

data . dat 

defines STAGS1 FOR002.dat file 
needed to run STAGS2 without 
having to rerun STAGS1 each time. 
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In the above statements, the file name 'data' is arbitrary and could be 
the Case Title. However, the file name extensions are not arbitrary and must 
be as defined above. 

2.2.4 SAMPLE STAGS RUN ON THE VAX 

Prior to execution, the user must have correctly prepared 
'data'.inp - data file for STAGS1 

and 

'data'. bin - data file for STAGS2 

user: RUNSTAGS <CR> 

system: Enter Root File Name. . . 

user: 'data' <CR> 

system: Print Log File? Y or N. 

user: Y for yes or N for no <CR> 

system: Run after 9PM? Y or N. 

user: Y for yes or N for no <CR> 

system: Is this a RESTART? Y or N. 

user: Y for restart or TP, N otherwise<CR>. 

system: Ok to submit STAGS run? Y or N. 

user: Y to submit or N to abort. <CR> 

STAGS2 execution on the VAX requires a series of new input records. The 
Thurston step/Transformation Processor (Ref. 3) referred to here is explained 
in Section 4. 

The total user input to the program comes in two parts: 

1. Input to STAGS 1 . This input block defines the model. It is 
completely unchanged from an ordinary model definition with the 
nonlinear static option turned on (INDIC-3). 

2. Input to STAGS2. This input file contains all necessary 
functions for a static, nonlinear analysis (including corotation 
and the Riks algorithm) . 
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The new input records TP-1, TP-2 and TP-3 contain the same data as the 
STAGS1 records A-l, C-l, and D-l, respectively. Data in the new records 
override the STAGS1 input records. This has the advantage that multiple 
solutions with differing load strategies can be run with the same STAGS1 
data. The present version of STAGS requires an additional record TP-4. 


TP-4 Secondary Path Parameter Record (see page 43 for parameter list) 

This record is used to select either the Thurston step/Transformation 
processor (TP) or the Riks algorithm, and when selecting TP, to specify the 
type of modal analysis desired and which mode to use for the shifted point. 
If NPATH is 0, STAGS2 will attempt to continue the solution along the current 
load path, just as would be done in an ordinary STAGS solution. If NPATH is 
1, a Thurston step will be attempted, and additional data from this record 
will be used. NEV specifies the number of bifurcation solutions to be 
calculated. NEQ specifies the number of modes to be used in the Thurston 
Equivalence Transformation. At present, NEQ cannot be greater than 1. NSOL 
describes the particular action of STAGS2 , accounting for multiple modes. 
For NSOL = 1 , it is assumed that only one mode is active, and a Thurston step 
is attempted using this mode for the shifted point. For NSOL=0 , the modal 
coefficient matrix is saved and the STAGS2 run is terminated. IE denotes the 
mode number to be selected in computing the shifted point. 


2.2.5 SAMPLE POSTP RUN ON THE VAX 

Prior to a run, the user must have correctly prepared 


1. 

'data'.inpp - data file for POSTP 


2. 

'data'.dat - data file created by 

a STAGS 1 run 

and 3 . 

'data'. sod - solution data file from the corre: 

user : 

RUNPOSTP <CR> 


system: 

Enter Root File Name... 


user : 

'data' <CR> 


system: 

Print Log File? Y or N. 


user : 

Y for yes or N for no <CR> 


system: 

Run after 9PM? Y or N. 


user : 

Y for yes or N for no <CR> 


system: 

Ok to submit POSTP run? Y or N. 


user : 

Y to submit or N to abort. 

<CR> 
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2.2.6 EXECUTING EAC ON THE MicroVAX 

To run EAC on the VAX, the user enters the command 
RUN EAC 

Entering this symbol initiates a command file which will submit an EAC 
analysis as a batch job. The user may directly access the executables for 
interactive execution: 

RUN DUAO : [ STAGS ] EAC . EXE 
and/or 

RUN DUAO : [ STAGS ] EAC PL . EXE 
The square brackets are needed! 

The user must use the following naming convention for the input files 
and resulting output files: 

defines the discrete FEM dk ta to be 
analyzed (TAPE43) 

defines the EAC execution data 7ile. 

(TAPE44) described as the ADDITIONAL 
INPUT FILE in Section 2. 
defines the intermediate file to which 
the generalized co-ordinates are written 
(TAPE45) . 

defines the data file for the plotting 
postprocessor EACPL and is created by 
EAC (TAPE46) 

defines the EACPL execution data file 
described in Section 2 (TAPE47). 
defines the EAC solution file, 
defines the log file for system 
responses . 


data. t43 
data. t44 

data . t45 

data. t46 

data . t47 

data. out 
data. log 
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Again, the file name 'data' is arbitrary . while the file name extensions are 
not arbitrary and must be defined as above. The corresponding naming 
convention on the Cyber is given in parenthesis. The file definitions are 
the same as for the Cyber version. 

2.2.7 SAMPLE EAC RUN ON THE MicroVAX 


Prior to the execution, the user must have valid data files as defined 
above . 


user : 

system: 

user : 

system: 

user: 

system: 

user : 

system: 

user: 

system: 


RUN EAC <CR> 

Enter Root File Name. 

'data' <CR> 

Print Log File? Y or N. 

Y for yes and N for no <CR> 
Run after 9 PM? Y or N. 

Y for yes and N for no <CR> 
Is this a RERUN? Y or N. 

Y for yes and N for no <CR> 
OK, Submitting EAC file. 


The results are returned in data. out and the graphics file as METAFILE.DAT. 
Currently the plots can be viewed only on TEK4109 terminals and copies made 
on hard copy machines. Alternatively, the plots can be printed (without 
viewing) on a printer such as the Post Script printer. To view the plots (on 
TEK4109 terminals), the user enters: 

SETDRV 409 <CR> 

M IRS HARE <CR> 

The system loads the metafile translator and displays a new prompt 'M>' 

M >S MF 1 METAFILE (assign the graphics file) 

M >D P 1 (to draw picture one) 

Press the spacebar to erase picture and continue. 

M >D P n (to draw the 'n'th picture) 

M >q (to terminate session) 

To make a hard copy, on TEK 4109 terminals, press the SCOPY button. 

To print the contour plots on the local Post Script printer in Building 1148, 
enter : 
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SET DRV PST <CR> 

MTRSHARE <CR> 

M >S MF 1 METAFILE 

M >D P FROM FIRST TO LAST (to draw all pictures) 

M >Q (to terminate session) 

This creates a data file called POSTS.DAT. The user may wish to rename this 
file. The command 'PSDI3000 fn' will submit this file as a batch job to the 
Post Script printer. 

NOTE: These instructions are only valid for the MicroVAX in Building 1148, 
Other users should refer to their systems coordinator. 

2.2.8 INTERACTIVE USAGE OF EAC ON THE MicroVAX 

Occasionally, the user may wish to execute EAC interactively on the VAX. 
This would certainly be necessary if the user wishes to evaluate the 
displacements , stress and moment resultants and transverse shear on a point 
by point basis. Prior to an interactive session, the user must assign the 
various data files used, to the corresponding I/O units. The following 
entries will make the necessary assignments: 

ASSIGN ' DATA' . T43 FOR043 
ASSIGN ' DATA' . T44 FOR044 
ASSIGN ' DATA' . T45 FOR045 
ASSIGN ' DATA' . T46 F0R046 
ASSIGN ' DATA' . T47 FOR047 
ASSIGN ' DATA ' . OUT FOR006 
To execute, enter 

RUN DUAO : [STAGS] EAC 

The results will be written to data. out. To generate the graphics file, enter 
RUN DUAO : [ STAGS ] EACPL 

Do not do this for single point evaluation sessions. 

2.2.9 MODIFYING EAC FOR LARGE GRIDS 

At the present time, the input and output data array dimensions are 
limited to a maximum of 50 rows and 50 columns each. To handle larger 
arrays, the user will have to access and modify a copy of the EAC source 
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code. This is available in the STAGS directory as a set of FORTRAN files. 
The command 

COPY DUAO: [STAGS. EAC . FOR] *.* *.* 
will copy the following files into the user's current directory: 

EAC . FOR 
NHM . FOR 
TRIMCT . FOR 
ERRANA. FOR 
EACS.OLB 
RUN EAC . COM 
MAKE_EXEC . COM 
EAC PLOT. FOR 

To modify the limits on the input (STAGS) grid, edit the FORTRAN files (files 
with FOR as the three letter extensions) and modify the values of SNX and SNY 
in the PARAMETER statements. SNX and SNY correspond to NRWS and NCLS 
respectively. To modify the limits on the output (user) grid, edit 
ERRANA. FOR and modify the values of NXU and NYU in the PARAMETER statements. 
In either case, EACPLOT.FOR must be modified. To make a new (local) 
executable , enter 
(2MAKE_EXEC 

This command file will compile and load the files to generate a local file 
called EAC. EXE. NOTE: Remember that the grid size cannot be made larger than 
what the computer can handle (memory limitations). 
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SECTION 3: THEORY 


This section contains the equations solved by the numerical analysis in 
EAC. The equations serve to define the problem actually solved by the 
analysis and to aid in comparing the output from EAC to the discrete finite- 
element results from STAGS. The main output is a set of generalized 
coordinates that allow computing a continuous solution for the plate 
equations. The continuous solution has the obvious advantage that stress 
resultants and deflections can be computed at any point on the plate and is 
not restricted to grid points in STAGS. However, the main advantage of the 
continuous solution is that a measure of its error can be computed. The 
output of the code includes this measure of the error in the continuous 
solution, which is also an indirect measure of the error in the finite 
element results in STAGS. 

In order to compute the continuous solution for a given plate from a 
shell structure, it is necessary to make some assumptions about the boundary 
conditions. In Ref. 4, the correction of the continuous solution by Newton's 
method is outlined along with removing the assumptions on the boundary 
conditions. The error measure computed in EAC is the starting point for this 
analysis, but correction of the results is beyond the scope of the analysis 
in the current version of EAC. 

3.1 CONSTITUTIVE RELATIONS 

The theory applies to orthotropic plates with the following constitutive 
relations : 
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where the strain-displacement relations for the imperfect plate are 
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3.2 EQUILIBRIUM EQUATIONS 


(3.3a) 

(3.3b) 

(3.3c) 


The equations that determine the transverse shear stress resultants are 
moment equilibrium equations 
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The in-plane equilibrium equations are 


(3.4a) 

(3.4b) 


N + N - 0 
x , x xy , y 

N + N - 0 
xy,x y,y 


The transverse force equilibrium equation is 


(3.5a) 

(3.5b) 


Q x,x + Vy ‘ N 3 (F ' w+ ^ } ” 0 (3.5c) 

The notation N^(F,W+W) is a shorthand notation 

N 3 (F,W+W) - N 3 (F,W) + N 3 (F,W) 

N,(F,W) * N W + 2N W + N W 
3 x , xx xy , xy y , yy 

3.3 DISPLACEMENT EQUATIONS AND BOUNDARY CONDITIONS 

The constitutive relations, strain-displacement relations, and moment 
equilibrium equations are substituted to reduce the three force equilibrium 


(3.6a) 

(3.6b) 


31 



equations, equations (3.5), to three partial differential equations in the 
three displacement components U,V, and W: 


L 11 (U) + L 12 (V) + N 1 (W ' W) + N 1 < w, * ) + ^(W.W) - 0 (3.7a) 

L 12 (U) + L 22 (V) + N 2 (W ’ W) + N 2 (W ' ; '' ) + N 2 (r ^’ W) = 0 (3.7b) 

L 33 (W) - N 3 (F,W) - N 3 (F,W) - 0 (3.7c) 


where 


L 11 (U) * A 11 U ,xx + A 66 U ,yy 
L 12 (f) - (A 12 + A 66 ) f jxy 


L 22 (V) 


A„ 0 V + A,, V 
22 ,yy 66 ,xx 


L--(W) -D,. W +2 (D 10 + 2 D,,) W + D 00 W 

33 11 ,xxxx 12 66 ,xxyy 22 ,yyyy 


(3.8a) 

(3.8b) 

(3.8c) 

(3.8d) 


N. (f,g) - A., f 


11 ,x 6 ,xx 


+ (A..„+A,,) f g + A f g (3.8e) 

12 66 ir 0 vi t hA v 0 \r\r 


,y 6 , xy 66 ,x & ,yy 


N 2 (f,g) - A 22 f (y g yy + <^ A i2 +A 66 ) f ,x S ,xy + A 66 f ,y s ,xx (3 ' 8f) 


In the finite - element analysis in the STAGS code, the plate shell unit 
is usually part of a larger structure. The shell unit must satisfy geometric 
and natural boundary conditions and/or continuity conditions between edges 
of the plate and the rest of the shell structure. In the EAC code, the STAGS 
results for displacements on the plate boundary become the discrete boundary 
conditions for a continuous solution of the nonlinear partial differential 
equations in the dependent variables U,V, and W, equations (3.7). 


3.4 CONTINUOUS SOLUTION FOR W 


The continous solution for W consists of two series solutions, W and 

h 

w with generalized coordinates computed from the discrete STAGS data. 

W - W_ + W_ (3.9a) 

1 h 

The functions W, and W' are chosen to satisfy two partial differential 
I E 

equations , 
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L 33 ( V * N 3 (F 0 >W 0 } + N 3 < F 0' W) 


(3.9b) 


L 33 (W e ) - 0 (3.9c) 

The zero subscript on the right side of equation (3.9b) denotes discrete 
STAGS data. The equation is only satisfied at the centroids of the 
rectangular elements as defined by the STAGS1 input on the number of rows and 
columns. The continuous solution for is computed by first using numerical 

harmonic analysis, Ref. 5, to compute the coefficients in a double sine 

series passing through the discrete values on the right-hand side of equation 
(3.9b), 

N 3 (F 0 ,W 0 +W) - 111 N 3mn sin m?rx/a sin mry/b (3.10a) 

The number of terms in the series is determined by M and N, the number of 

rows and columns in the STAGS1 model ). The generalized coordinates in the 

continuous solution W T are the Fourier coefficients A in the series 

I mn 


W 


I 


s-f 

m-1 


l-l 

i-l 


A sin nwrx/a 
mn 


sin n;ry/b 


(3.10b) 


Substituting equations (3.10) into equation (3.9b) determines the numerical 
values of the A 

mn 



N 


3mn 


[D n m 4 + 2 <D 12 + 2 D 66 ) m 2 nV + D^nV] 


(3.10c) 


In general, the continuous solution W does not satisfy the plate boundary 
and/or continuity conditions. The solution is added to to obtain a 
better fit of the discrete results on the boundary. The solution is 
assumed as a series that satisfies equation (3.9c) term by term. 
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(3.11) 


“e ' «oo * |-i IV (,) cos p? * V ,) sin Pfl 


+ 2 - [W (£) cos qrj + W (£) sin qr; ] 
q-1 qc qs n 


The first term is a polynomial with undetermined coefficients 


W 00 - a 00 + a 10 x + a 01 y + a 20 x2 + a 02 y2 


(3.12) 


Each of the functions W (rj) , W (r/) , W (£) . and W (£) satisfies a fourth 

pc ps qc qs 

order, homogeneous, ordinary differential equation and, consequently, 

contains four constants of integration that make up the set of generalized 

coordinates for W_. 

E 


3.5 BOUNDARY CONDITIONS FOR W 


The undetermined coefficients in and the constants of integration in 

the functions W (rj) , V (rj) , W (?)» and W (O are computed from discrete 
pc ps qc qs 

STAGS output on the plate boundary. The output for W, W on sides with 

i x 

constant x, and W on sides with constant y at nodal points become the 

t y 

discrete boundary conditions for the continuous solution for W. 

Accepting the STAGS results as correct on the plate boundary is an 
assumption in deriving the continous solution in EAC . After all the 

generalized coordinates in the continuous solution for U , V , and W are 
computed, the effects of the assumption can be examined. However, correction 
of errors on the boundary of the plate is beyond the scope of the theory 
programmed in EAC. 


The constants of integration in the product solutions in W are computed 

E 

by a least- squares analysis of the equations, 
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X “0 and x=-a 


( 3 . 13a t b) 


u - u . u 
E 0 I 

“ w n ' W T 

E ,x 0 , x I , x 

u * u - u 
E 0 I 


W_ - W A W_ 

E.y O f y I ,y 


x -0 and x*a 
y =0 and y=b 
y =0 and y=b 


( 3 . 13c , d) 
(3 . 14a,b) 
( 3 . 14c , d) 


For example, four equations containing the four constants of integration 
appearing in the function W^^ry), where r is an integer, are generated by 

multiplying each of the equations (3.13) by cos rr? and integrating from 0 to 

2n . Terras in W ( 17 ) for q^r drop out of the four integrals along with all 
qc 

the functions W (n) . However, the integrals determine nonzero coefficients 
qs 

multiplying the constants of integration appearing in the functions W (tj) 

pc 

and W (ry) . 
ps 


In a corresponding fashion, four integrals are written for each 
trigonometric function cos qr? , sin qrj , cos pf , and sin p£ appearing in W^,. 

Integrals are also written containing nonzero multipliers of the undetermined 
coefficients in W qq . Four of these latter integrals are integrals of 

equations (3.13c), (3.13d), (3.14c) and (3.14d). An additional equation 

containing the coefficient a^ is derived by summing the four integrals of 

equations (3.13a), (3.13b), (3.14a) and (3.14b). 

The integrals derived from the left-hand sides of equations (3.13) and 
equations (3.14) are equated to corresponding integrals of the right hand 
sides. Since the notation refers to discrete STAGS data at nodal points, 

the integrals containing and its partial derivatives are replaced by 

numerical quadrature formulas. 

In order for the numerical quadrature formulas to be meaningful, there 
is an upper limit p-P and q-Q for the number of terms in the series for W^,, 

equation (3.11). The upper limits P and Q must satisfy the inequalities (see 
ref. 5) 
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2P + 1 < M 


(3.15a) 


2Q + 1 < N (3.15b) 

where M and N are the number of rows and columns in the STAGS data. The 
number of terms, P and Q, correspond to the FORTRAN variables, MP and NQ on 
card EAC-3 in the file TAPE44, and M and N are NROWS and NCOLS in the input 
file for STAGS 1 . 

The integrations defined above result in a set of linear algebraic 
equations for a least- squares fit of the STAGS boundary data on W. The 
equations are solved to compute numerical values for (8P+8Q) constants of 
integration and for five undetermined coefficients in the polynomial W qq ■ In 

the computations that follow in EAC , the continuous function i s used 

wherever W or its partial derivatives appear in the equations. 

3.6 CONTINUOUS SOLUTIONS FOR U AND V 

The continuous solutions for U and V are computed following the same 
general plan as the solution for W. A double Fourier series is used to 
generate solutions and of the in-plane equilibrium equations, equations 

(3.7a,b). Solutions of 


11< U E> 

+ L 12 (V E> ' ° 

(3.16a) 

12< U E> 

+ L 22<V ‘ 0 

(3.16b) 


are added so that continuous solutions of U and V are of the form 


U " U ! + U E 

(3.17a) 

V - V ! + V E 

(3.17b) 


Constants of integration appearing in the solutions U £ and V^, are 

computed from a least- squares fit of discrete STAGS data for values of U and 
V at node points. 
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The specific forms for U and V are 

U T - 2 A 2 n A- cos m£ cos nr; + 2 - 2 -.B- sin mf sin nr? 

I m— 0 n-0 lmn s ' m-1 n-1 lmn , s ' 

+ 2 A 2 - C- cos mf sin nr? + 2 - 2 A D.. sin mf cos nr? (3.18a) 

m=0 n-1 lmn m-1 n-0 lmn 

V - 2 A 2 A A 0 cos m£ cos nr? + 2 - 2 sin m£ sin nr? 

I m-0 n-0 2mn m-1 n-1 2mn s 

+ 2 A 2 - C 0 cos mf sin nr? + 2 - 2 A D 0 sin m£ cos nr? (3.18b) 

m=0 n-1 2mn m-1 n-0 2mn s ' 

U E * a l + b l X + C l y + p-l^ F lp^ C0S + G lp (r?) Sin 

+ q**l [ H l q (^ ) cos qr? + I lq^ ) sin qr} J (3.19a) 

V E ~ a 2 + b 2 X + C 2 y + p-l^ F 2p (r?) C0S P ^ + G 2p (??) Sin p ^ 

+ £^[^(0 cos qr? + I 2q^ ) sin qr ^ (3.19b) 


For solution of equations (3.16), the functions F . (r?) and F 0 (r?) satisfy the 

lp 2p 

pair of ordinary differential equations 


A 66 ^ F lp ' p2 A 11 F lp + (A 12 + A 66 
' (A 12 + A 66 ) F lp + A 22 ^ F 2p ‘ p2 


> ^ P F 2p 
A 66 F 2p 


0 

0 


(3.20a) 

(3.20b) 


where primes denote differentiation with respect to r? . The subscripted 
functions G, H, and I satisfy similar coupled ordinary differential 
equations. The complete solution of equations (3.20) contains four constants 
of integration. The constants of integration are determined by boundary 
conditions on U and V. 


V. 



U. 

V. 

u. 

V 


x-0 and x-a 
x=0 and x-a 
y-0 and y-b 
y=0 and y-b 


( 3 . 21a ,b) 
(3 . 21c , d) 
( 3 . 22a , b) 
( 3 . 22c , d) 
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The solution for the undetermined coefficients and constants of integration 
in U_ and V is analogous to the solution for W In particular, the number 

of terms in the series solutions must satisfy the inequalities in equations 

(3.15). The FORTRAN variables MT and NT on card EAC-3 in the file TAPE44 

correspond to the number of p and q terms in U and V The upper limits on 

Eh 

p and q must again satisfy the inequalities of equations 3.15. 

3.7 OUTPUT AND ERROR ANALYSIS 

The leas t- squares solution for the generalized coordinates in and 

completes the analysis to determine a continous solution for the plate 

displacements U, V, and W. The stress resultants N , N , N , M , M , M , 

r x y xy x y xy 

Q x , and are then computed from the displacements and their derivatives. 

The output from the continuous solution can be computed at any given point or 
over a rectangular grid (see EAC cards in NPOSTPX) . 

A measure of the error in the continuous solution derived from the 

finite element results is E^(x,y), the residual error in satisfying the 

transverse equilibrium equation, equation (3.7c). 


E 3 * L 33 (W) ' N 3 (F,W+W) 


(3.23) 


where W is the continous solution defined by equations (3.9) and the in-plane 
stress resultants appearing in the bilinear operator are computed from the 
continous solution. The output from EAC normalizes E^ by dividing the 

residual error by the maximum absolute value computed for the linear operator 


Away from the plate boundary, the error is a good measure of the 
accuracy of the results from the continuous solution and represents a result 
that is not available in a f ini te - element solution based on direct 
minimization of a variational principle. Near the plate boundary, the error 


38 



contains a 


"Gibb's phenomenon" because 
and the imperfection W. On the boundary 

l 33 (w i ) - 0 

and 

E 3 - - N 3 (F,w E ) 


while away from the boundary the general 


of using a double sine series for W^ 


error, from equations (3.9), is 


E 3 * N 3 (F 0’ W 0 + ^ } ’ N 3< F ’ W+ ^) (3.24) 

On the plate boundary, other measures exist for examining the accuracy 
of the continuous solution. On boundaries, where STAGS1 specifies constant 
displacements, the output from EAC can be compared directly to the STAGS 
output at nodal points. When the maximum number of linearly independent 

generalized coordinates defined in equations (3.15) is used in the continuous 
solution, the least- squares results for displacements will agree with the 
STAGS results at nodal points on the boundary except at the four corners of 
the plate. The continuous solution can also be checked against displacement 
boundary conditions between nodes. 

In addition to displacements, stress resultants are also computed on the 
boundary in EAC. This allows checking the error in the stress resultants. 
For example, on a s imply- supported edge perpendicular to the x-direction, the 
bending moment resultant, M^, should vanish. The analysis in EAC matches the 

rotation, W , on this boundary and the residual M p on a simply-supported 
, x x 

edge is a measure of error although it will not be clear whether the error is 

in the continous solution or in the edge rotations from the finite-element 

solution. 

When continuity relations are specified for the edge of a plate at a 
juncture with other shell units, the accuracy of the results from STAGS and 
EAC are more difficult to assess. The STAGS results will satisfy continuity 
of displacements at the common nodes of the shell units at a juncture. The 
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continuous solution can also be used to check continuity between nodes along 
a common boundary. 

When a number of shell units make up a statically indeterminate 
structure, it is of interest to determine the load paths in the structure. 


EAC computes output for in-plane 
x-direction 



and the y-direction, 

F - - [ N dy 

yi J o *y 

F - f N dx 

y2 J 0 y 

F - [ N dy 

y3 J 0 x y 

F , - - f N dx 

y 4 J 0 y 


force resultants on the edges, in the 


x-0 

(3.25a) 

y-b 

(3.25b) 

x-a 

(3.25c) 

y-0 

(3 . 25d) 


x-0 

(3.26a) 

i 

y-b 

(3.26b) 

x-a 

(3.26c) 

y-0 

(3 . 26d) 
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SECTION 4: TP USER'S INSTRUCTIONS 


A new processor called TP (transformation processor) is available. TP 
allows the user to obtain a nonlinear solution on an equilibrium path 
different from the primary equilibrium path. This processor examines the 
solution on the primary path and for a fixed set of input parameters will 
compute solutions on another path, well into the postbuckling range. Ref. 
(2) outlines the numerical analysis for the TP processor and contains two 
example cases. 

The Bifurcation Point Processor, TP, is an executable program which 
replaces STAGS2 in the normal execution sequence (STAGS1, STAGS2) . TP can 
perform any of the following steps in one execution: 

1. Nonlinear collapse analysis from any specified restart solution to a 
specified maximum load level. The restart solution may be on a primary or a 
secondary path. The Riks path length method is used here. This step is the 
same as a normal STAGS2 execution using the Riks method. 

2. Nonlinear collapse analysis as in (1) above with the additional 
computation of the closest bifurcation modes whenever a matrix factorization 
is performed. The bifurcation modes are saved on the restart file TAPE22 and 
are thus available in later execution steps as initial estimates or for use 
in step 3. 

3. Secondary solution computations using Thurston's method. This type of 
execution can be performed as a separate step or can sometimes be preceded 
with computations described in (2) above. 

TP requires a short input file to control execution of the different 
options available to the user. The contents of this data file and the 
execution of TP differ depending on the computer, Cyber or MicroVAX. 

4.1 WORKING ON THE CYBER 

The input data are read in free format as described for the STAGS1 input 

file . 
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TP- 1 Case Title 


The case title may contain any text up to 72 characters. 

Variable Description 

TITLE Case title 


TP-2 Load Multiplier Record. 

This record is the same as STAGSC-1 record C-l. 


Variable 

STLD1 

STEP1 

FACM1 

STLD2 

STEP2 

FACM2 

ITEMP 


INERT 


Description 

Starting load factor for system A. 

Load factor increment for system A. 
Maximum load factor for system A. 

Starting load factor for system B. 

Load factor increment for system B. 
Maximum load factor for system B. 

0 - No thermal loads . 

1 - Temperatures referred to system A. 

2 - Temperatures referred to system B. 

0 - No gravitational/centrifugal effects. 

1 - Grav/cent loads in system A. 

2 - Grav/cent loads in system B. 


TP- 3 Strategey Parameter Record. 

This record is the same as the STAGSC-1 record D-l. 


Variable 
I START 

NSEC 

NCUT 

NEWT 


NSTRAT 


Description 

0 - Begin new case . 

N - Restart from Nth load step. 

No. of CPU seconds of execution (approx). 
Total no. of times the step size may be cut. 
Total no. of refactorings allowed. 

NEWT < 0 Refactor every NEWT load steps. 
NEWT**- 20 True Newton method. 

0 - Initial estimate from extrapolation. 

1 - Previous soln. used as initial estimate. 
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-1 - Path length used as independent param. 

DELX Error tolerance . 

(if DELEX-0, default DELEX-.001 is used) 

WUND Relaxation factor. 

(if WUND=0, the program chooses WUND) 

TP-4 Secondary Path Parameter Record . 

The parameters in this record control the branching from the primary path to 
a secondary path. After solutions on a secondary path have been obtained and 
saved on the restart file, the analysis may be continued with the standard 
"primary" path method using a secondary path solution as the restart load 
step . 


Variable 

NPATH 

NEV 

NEQ 


NSOL 


IE 


Description 

0- Primary path (Riks' method). 

1- Secondary path (Thurston's method). 

0- No bifurcation solutions will be found. 

1- Find N bifurcation modes. 

0-Do not continue to bifurcation path. 

N- Equivalence transformation using N modes. 
(NOTE: NEQ cannot exceed NEV) 

0- Save reduced system on TAPE31; stop 

1- Save reduced system on TAPE31; solve for 
load factor PA with QFAK=0 using 
equivalence transformation. 

-1-Read PA, QFAK (TP-6); obtain solution using 
equivalence transformation . 

Number of bifurcation mode used to compute 

shifted point. (Not used if NPATH or NEQ = 0.) 


TP- 5 Secondary Solution Path Reduced System Record. 

This record is included only if NPATH=1 and NEQ>0 on the TP -4 record. STEP 
specifies the amplitude of the eigenvector of the tangent stiffness matrix 
chosen to produce a secondary solution. (The STAGS output labels these 
eigenvectors as "buckling modes"). The "first shifted point" is obtained by 
adding STEP times the IE'th mode to the last solution on the primary path. 
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The eigenvectors have been normalized so that the largest component is 1. 
DELTA determines the perturbations of the NEQ eigenvectors from the "first 
shifted point". A reduced polynomial system is thus obtained to provide for 
numerical solutions to the complete nonlinear problem. 


Variable 

STEP 


DELTA 


Description 

Amplitude of bifurcation mode IE used to 
construct the shifted point 
(i.e. , XO - X+STEP*EV (IE)) . 

Coefficient of perturbations used for 
reduced system. 


TP- 6 Eigenvector Amplitude Record. 

This record is included only if NS0L=-1 on the TP-4 record. The equivalence 
transformation determined by NEQ eigenvectors is used with the load factors 
PA and PB and modal amplitudes QFAK to obtain a solution in the full discrete 
system. 


Variable 

PA 

PB 

QFAX(I) 
(1=1, NEQ) 


Description 

Load factor for system A for secondary solution. 
Load factor for system B for secondary solution. 
Amplitudes of eigenvectors for 
secondary solution. 
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4.1.1 STANDARD TRANSFORMATION PROCESSING (TP) 

JOB. . . 

USER. . . 

CHARGE. . . 

DELIVER. . . 

FILESET , IA-CPF , DA-LPF . 

GET , STAGS/UN-STAGS , ST-CPF . 

GET , TAPE20*pfn . 

BEGIN(TP, STAGS) 

STAGS 1. 

STAGS 2. 

STAPL. 

PLOT . VARIAN , FF 
E-O-R 

STAGS1 data 
E-O-R 

TP data 
E-O-R 

STAPL data 
E-O-F 


If user-written subroutines are included, it will be necessary to reload 
the STAGSC-1 processors using the procedure RLTP . This process is similar to 
that described in section 2.2.1 for RLBIG and RLPOST. 
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4.2 WORKING ON THE MicroVAX 


On the MicroVAX, TP is still called STAGS2 and always requires input. 
The input to TP was described in Section 2 under the heading SAMPLE STAGS RUN 
ON THE VAX. To select and use the TP algorithm, three of the input cards are 
modified. 

On TP-2 STEP is replaced by the value of the estimated modal amplitude. 

On TP-3 NSTRAT must be set to -1 for the TP option to be used. 

On TP-4 select NPATH to be 1. If NPATH - 1, a Thurston step will be 
attempted, and additional data from this record will be used. NEV specifies 
the number of eigenvectors to calculate. NEQ specifies the number of modes 
to be used in the Thurston Equivalence Transformation. NSOL must be selected 
as unity. 



APPENDIX A: NOMENCLATURE 


Analysis 

Variable 

FORTRAN 

Name 

Description 

A 11’ A 12 ,A 22’ A 66 

All , A12 , A22 , A66 

plate extensional stifnesses 

A , N, 
mn 3mn 


Fourier sine-sine series 
coefficients 

a,b 

A, B 

dimensions of rectangular plate 

D ir D 12’ D 22’ D 66 

Dll , D12 , D22 , D66 

plate bending stiffnesses 

e , e , e 
x y xy 


in-plane strains 

E 

E3 

residual error 



dummy variables in linear 
operators 

L tj ( ) i-1,2; j-1,2 


linear operators 

M , M , M 
x y xy 

MX,MY,MXY 

moment stress resultants 

N f N , N 
x y xy 

NX , NY , NXY 

in-plane stress resultants 

x,y 

XX(I) , YY( J) 

Cartesian coordinates 

u,v,w 

u t v,w 

plate displacements 

w 

WI 

initial out-of-plane 
imperfection 

0=a/b 

BETA 

plate aspect ratio 

£**2nx/a, rj*-2ny/b 


dimensionless coordinates 
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APPENDIX B: EXAMPLE 


A stiffened plate assembly (Ref. 6) is shown in Fig. 1. The general 
purpose finite element code STAGS is used to analyse this plate assembly in 
the postbuckled state. The plate assembly is simply supported on two edges 
and symmetric on the other two edges. A compressive load twice the critical 
load is applied on the simply supported edges. Data for this sample problem 
are given in Table 1. Assuming the finite element results are correct on the 
edges of the shell unit, the error analysis code EAC is used to obtain 
numerical results for shell unit 3, a plate that is bounded on two sides by 
a juncture with a stiffener and another plate shell unit. Sample data files 
are given here for the STAGS and the EAC runs on the Cyber and on the 
MicroVAX. Sample contour plots generated by EACPL are also included in 

Figs. 2-13. 

Note: Currently EAC does not run with the corotation option. 

WORKING ON THE CYBER 

For a line by line explanation of the STAGS input data cards, the user 
is referred to the STAGS users manual (Ref. 1). These data cards indicate 
that STAGS is to use the non-corotational procedure and nonlinear analysis. 
There is a single loading system compressing the panel and analysis is to be 
done from a starting load factor of 0.4 to a final load factor of 1.4 in 

steps of 0.2. There are 5 shell units. The connections between the shell 

units are given on the G-l cards and the partial compatibility constraints 

for uniform end shortening are given on the G-2 cards. Material and wall 

data are given on the I and K cards. Shell unit details are given on the M 
cards . 

STAGS 1 INPUT DATA 

5 BRANCH STIFFENED PANEL MODEL 
3110001001 $B-1 

5 0 0 4 10 0 0 $B-2 

102 $B- 3 
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0.4 0.2 1.4 $C-1 

0 5000 0 -20 0 1.0E-4 $D-1 

17,5 17,3 17,7 17,3 17,5 $F-1 

C 

C DEFINE CONNECTIONS BETWEEN SHELL UNITS 

C 

1234 $G-1 

1224 $G- 1 

3244 $G- 1 

3254 $G-1 

C 

C DEFINE PARTIAL COMPATIBILITY CONSTRAINTS 

C FOR UNIFORM END SHORTENING 

C 

3131 1101 $G-2 

3131 2101 $G-2 

3131 3101 $G-2 

3131 4101 $G-2 

3131 5101 $G-2 

3 17 3 1 1 17 0 1 $G-2 

3 17 3 1 2 17 0 1 $G-2 

3 17 3 1 3 17 0 1 $G-2 

3 17 3 1 4 17 0 1 $G-2 

3 17 3 1 5 17 0 1 $G-2 

C 

C MATERIAL AND WALL DATA 

C 

1 $ 1-1 

1.32E7 0.3 4.54E6 0.101 0.0 9.8E6 0.0 $1-2 

111 $K- 1 

1 0.084 $K-2 

211 $K- 1 

1 0.058 $K-2 

C 

C UNIT 1 
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c 

2 3 $M-1 

0.0 15.0 0.0 2.5 
0.0 0.0 0.0 $M-4 

0.0 2.5 0.0 $M-4 

15.0 2.5 0.0 $M-< 

1 $M-5 

411 $ 

1614 $P-1 

0000 $Q-1 

1 5 $R-1 

C 

C UNIT 2 
C 

23 $M- 1 

0. 15.0 0.0 1.352 
0.0 2.5 0.0 $M-4 

0.0 2.5 1.352 $M 
15.0 2.5 1.352 $1 

2 $M- 5 

411 $ 

1616 $P-1 

0000 $Q- 1 

1 5 $R-1 

C 

C UNIT 3 
C 

2 3 $M- 1 

0.0 15.0 0.0 5.0 
0.0 2.5 0.0 $M-4 

0.0 7.5 0.0 $M-4 

15.0 7.5 0.0 
1 $M- 5 

411 $ 

1616 $P-1 

1000 $Q- 1 


$M-2A 


$M-2A 

4 

:-4 


$M-2A 
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1 3 $Q-2 

8000.0 11130 $Q-3 

-8000.0 1 1 17 3 0 $Q- 3 

0.0-11930 $Q- 3 ELIMINATE U RIGID BODY MOTION 

1 5 $R-1 

C 

C UNIT A 
C 

2 3 $M- 1 

0.0 15.0 0.0 1.352 $M-2A 

0.0 7.5 0.0 $M-4 

0.0 7.5 1.352 $M-4 

15.0 7.5 1.352 $M-4 

2 $M-5 

411 $ 

1616 $P-1 

0000 $Q- 1 

1 5 $R- 1 

C 

C UNIT 5 
C 

2 3 $M- 1 

0.0 15.0 0.0 2.5 $M-2A 

0.0 7.5 0.0 $M-4 

0.0 10.0 0.0 $M-4 

15.0 10.0 0.0 $M-4 
1 $M- 5 

411 $ 

1416 $P-1 

0000 $Q- 1 

1 5 $R- 1 

TP DATA FILE 


Whereas STAGS2 does not normally read input data, TP reads a short input 
file to control the execution of the different options available to the user. 
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This input file also permits the user to change certain data records which 
may have been read in STAGS1. On TP-2, the initial and final load factors 
are specified as the load factor at which TP is to be attempted. The data 
below are representative. 

5 BRANCH STIFFENED PANEL MODEL. TP RUN $TP-1 
1.0 0.1 1.0 $TP-2 

4 8000 2 -20 1 0 $TP- 3 

12111 $TP-4 

0.04 0.1 $TP- 5 

DATA FOR RESTART FROM TP 

After solutions on a secondary path have been obtained and saved on the 
restart file, the analysis may be continued with the standard "primary" path 
method using a secondary path solution as the restart load step. In the 
STAGS1 data file, modify the C-l and the D-l records. 

1.204 0.2 2.0 $C-1 

5 8000 0 -20 1 1.0E-4 $D-1 

WORKING ON THE MicroVAX 


STAGS 1 DATA 

Same as for Cyber. 

STAGS2 DATA 

There is no corresponding data file on the Cyber. The first three input 
records TP-1, TP-2 and TP-3 contain the same data as the STAGS1 A-l, C-l, and 
D-l records, respectively. Data in the new records override the STAGS1 input 
records. 

5 BRANCH STIFFENED PANEL MODEL $TP-1 
1.0 0.1 1.0 $TP-2 

0 8000 880 1.0E-4 $TP-3 

1 $TP-4 
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TP DATA 


When the Thurston Equivalence Transformation is to be used, the data 
item STEP on record TP-2, now takes the value ot the modaL amplitude. The 
value of NSTRAT on TP- 3 is -1. The value of NPATH on TP-4 is 1 and 

additional data from this record will be used. In the VAX version, STAGS 
will continue the analysis along the secondary path till a desired load 
factor (as specified in item 3 of TP-2 card) is achieved. On the Cyber, the 
user has to restart from the load factor to which TP converges. 

3 BRANCH STIFFENED PANEL MODEL $TP-1 
1.0 0.04 2.0 $TP- 2 
1 8000 2 -20 -1 0 $TP- 3 
12112 $TP-4 

The modal amplitude has been specified to be 0.04; the equivalence 
transformation uses 2 vectors with eigenvector 2 as the initial direction. 
Two eigenvectors (denoted as "buckling modes" in the STAGS output) will be 
computed. Start analysis from step 1 of the previous run. Continue the 
analysis along the selected path to a load factor of 2.0. 

RESTARTS 

A STAGSC-1 analysis can be restarted from any load step for which the 
solution has been saved on a restart file during a previous execution. 

Example : 

5 BRANCH STIFFENED PANEL MODEL $TP-1 
1.220 0.4 2.622 $TP-2 

5 8000 0 -20 1 0.05 $T?-3 
1 $TP -4 

The analysis is restarted from the 3th load step with a starting load factor 
of 1.22 and a final load factor of 2.622, with a load step of 0.4. The 
execution time is expected to be less than 8000 seconds. 
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EAC INPUT DATA 


EAC is used to obtain numerical results for plate component 3 that is 
bounded on two sides by a juncture with a stiffener and another plate 
component. The input data indicate that this is a new case, that results are 
to be generated over a grid defined in the next record using 8 terms in X and 
3 terms in Y for both the inplane and the transverse equations' homogeneous 
solutions. Printout of input data and output results is requested, and the 
analysis is for a plate with no defined initial geometric imperfections. 

SHELL UNIT 3 OF 5 ELEMENT STIFFENED PANEL 

0 $EAC-2 

. 126176E+07 . 378527E+06 .936758E+06 .381360E+06 $EAC-2A 
. 741912E+03 . 222574E+03 .550814E+03 .224240E+03 $EAC-2B 
08383 $EAC-3 
33 13 $EAC-4 

1 0 $EAC-5 
0 $EAC-6 

EACPL INPUT DATA 


The following are the input data for the contour plotting postprocessor. 
The input grid is 17 by 7, the output grid is 33 by 13, the contour plots are 
to be scaled to 40% of the plate dimensions, and the contour limits are user 
specified in 12 EACPL- 2 records. Note: When contour data are user input, all 
12 EACPL-2 records must be present. 

17 7 33 13 0.4 1 $EACPL- 1 
-0.01 0.01 0.001 $ EACPL- 2 

-0.001 0.001 0.0001 $ EACPL- 2 

-0.1 0.1 0.01 $ EACPL- 2 

-2000.0 0.0 50.0 $EACPL- 2 
-500.0 500.0 50.0 $EACPL-2 
-50.0 50.0 5.0 $ EACPL- 2 

-30.0 30.0 2.0 $ EACPL- 2 

-30.0 30.0 2.0 $ EACPL- 2 
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-30.0 30.0 2.0 $EACPL-2 
-80.0 76.0 4.0 $EACPL- 2 
-80.0 76.0 4.0 $EACPL-2 
-0.02 0.019 0.001 $EACPL- 2 
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Table 1. Data for sample problem. 


Property 


Component 


Skin Stiffener 


A lx , kips/in 

12^2 

871 

A 12 , kips/in 

379 

261 

A 22> kips/in 

937 

647 

A 66> kips/in 

381 

263 

D 1L , in. -lb 

742 

244 

D 12 , in. -lb 

223 

733 

D 22> in. -lb 

551 

181 

D 66 . in. -lb 

224 

738 

Thickness, in. 0 

.084 

0.058 

Panel length, a - 

15.0 in. 


Middle plate width 

- 5.0 in. 


Outer plate widths 

~ 2.5 in. 


Stiffener height 

- 1.4 in. 



Load, P- 15,792 lb. 
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(c) Finite element grid. 


(d) Postbuckled deflection pattern. 


Fig. J . Exaasle of a plate caeponent from a 
stiffened panel. 
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(c) Continuous Least Squares Approximation 
Fig. 3 . In-plane Deflection, V, Contour Plots For Plate 2. 
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(b) Continuous Approximation (Ref. 4) 



(e) Continuous Least Squares Approximation 
Fig. 4. Transverse Deflection, V, Contour Plots For Plate 2. 
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(a) Discrete FEM Solution 











ut.i r'l'T/-! rr.T, I 


(b) Continuous Approximation (Ref. 4) 





i .rr i 

I .L-4-H'* 

1 ! 1 1 
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(c) Continuous least, Squares Approximation 
'ess Resultant, , Contour Plots For ruate 2. 
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(a) Discrete FEM Solution 



(b) Continuous Approximation (Ref. 4 ) 



(c) Continuous Least Squares Approximation 
Fig. 6. Stress Resultant, N , Contour Plots for Plate 2. 
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(a) Discrete FEM Solution 



'’ig. 7. Stress Resultant, » Contour Plots for Plate 2. 


ORfGiim PAGE ,'S 

OF POOR QUALITY 
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(a) Discrete F2H Solution 



(b) Continuous Approximation (Ref. 4 , 



(c) Continuous Least Squares Approximation 
Pig. S. Moment Resultant, M. , Contour Plots for Plate 2. 
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(b) Continuous Approximation (Ref. 4) 



(0) Continuous Lea so Squares Approximation 
Fig. 9 .. Moment Resuioant, M , Contour Plots for Plate 2. 
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(a) Discrete FEM Solution 



(t) Continuous Approximation (Ref. 4) 
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(a) Continuous Approximation (Ref. 4) 



(b) Continuous Least Squares Approximation 
rig. 11. Transverse Shear, Q , Contour Plots for Plate 2. 
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(a) Continuous Approximation (Ref. 4 ) 



(b) Continuous Least Squares Approximation 
Fig. 12. Transverse Shear, Q , Contour Plots for Plate 2. 


69 





I 

i 

ORIGINAL FI'3’1 «3 

OF. POOR QU/.UTY 



Fig. 13. Estimate Error in Transverse Equilibrium Equation. 
(Error =(L jf ^)/^ 33max ) 
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